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Abstract. Massive stars played a key role in the early evolution of the Universe. They formed 
with the first halos and started the re-ionisation. It is therefore very important to understand 
their evolution. In this paper, we describe the strong impact of rotation induced mixing and 
mass loss at very low Z. The strong mixing leads to a significant production of primary ^*N, 
^•^C and ■^^Ne. Mass loss during the red supergiant stage allows the production of Wolf-Rayet 
stars, type Ib,c supernovae and possibly gamma-ray bursts (GRBs) down to almost Z = Q for 
stars more massive than 60 Mq. Galactic chemical evolution models calculated with models of 
rotating stars better reproduce the early evolution of N/O, C/0 and ^^C/^'^C. We calculated 
the weak s-process production induced by the primary ^^Ne and obtain overproduction factors 
(relative to the initial composition, Z — 10~^) between 100-1000 in the mass range 60-90. 

Keywords. Stars: mass loss. Population II, rotation, supernovae, Wolf-Rayet, Galaxy: evolu- 
tion, gamma rays: bursts 



1. Introduction 

Massive stars {M > 10 Mq) started forming about 400 millions years after the Big 
Bang and ended the dark ages by re-ionising the Universe. They therefore played a 
key role in the early evolution of the Universe and it is important to understand the 
properties and the evolution of the first stellar generations to determine the feedback 
they had on the formation of the first cosmic structures. It is unfortunately not possible 
to observe the first massive stars because they died a long time ago but their chemical 
signature can be observed in low mass halo stars (called EMP stars), which are so old 
and metal poor that the interstellar medium out of which these halo stars formed are 
thought to have been enriched by one or a few generations of massive stars. Since the 
re-ionisation, massive stars have continuously injected kinetic energy (via various types 
of supernovae) and newly produced chemical elements (by both hydrostatic and explosive 
burning and s and r processes) into the interstellar medium of their host galaxies. They 
are thus important players for the chemo-dynamical evolution of galaxies. Most massive 
stars leave a remnant at their death, either a neutron star or a black hole, which often 
lead to the formation of pulsars or X-ray binaries. 

The evolution of stars is governed by three main parameters, which are the initial 
mass, metallicity {Z) and rotation rate. The evolution is also influenced by the presence 
of magnetic fields and of a close binary companion. For massive stars around solar metal- 
licity, mass loss plays a crucial role, in some cases removing more than half of the initial 
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mass. Internal mixing, induced mainly by convection and rotation also significantly af- 
fect the evolution of stars. The properties of non-rotating low-Z stars are summarised in 
Sect. 2 of Hirschi et al. (2008). In short, at low Z, stars are more compact, usually have 
weaker winds (see however IPulstilnik et al. 2008P and are more massive (Bromm & Loeb 
2003, [Schneider et al. 2006|) . The fate of non-rotating massive single stars at low Z is 



summarised in 



Heger et al. (2003) and several groups have calculated the corresponding 



stellar yields ( |Heger fc Woosley 2002l|Chieffi fc Limongi 20041 [Tominaga et al. 200"7| . In 
this paper, after discussing the impact of rotation induced mixing and mass loss at low 
we present the implications for the nucleosynthesis and for galactic chemical evolution 
in the context of extremely metal poor stars. 



2. Rotation, internal mixing and mass loss 

Massive star models including the effects of both mass loss and especially rotation 
better reproduce many observables around solar Z . For example, models with rotation 
allow chemical surface enrichments already on the main sequence (MS), whereas with- 
out the inclusion of rotation, self-enrichments are only possible during the RSG stage 
( Heger fc Langer 2000[[Meynet fc Maeder 2000 ). Rotating star models also better repro- 
duce the WR to O type star ratio and also the ratio of type Ib-|-Ic to type II supernova 
as a function of metallicity compared to non-rotating models, which underestimate these 
ratios (see contribution by Georgy in this volume and [Meynet fc Maeder 2005D . The 
models at very low Z presented here use the same physical ingredients as the successful 
solar Z models. The value of 300 kms~^ used as the initial rotation velocity at solar 
metallicity corresponds to an average velocity of about 220 km s^^ on the main sequence 
(MS) which is close to the average observed value. One of the first surveys of OB stars 
(and in particular their surface velocity) was obtained by Fukuda (1982) and the most 
recent surveys are listed in Meynet et al. (2008). It is unfortunately not possible to ob- 
serve very low Z massive stars and measure their rotational velocity since they all died 
a long time ago. Higher observed ratio of Be to B stars in the Magellanic clouds com- 
pared to our Galaxy ([Maeder et al. 19991 |Martayan et al. 2007[ ) could point out to the 
fact the stars rotate faster at lower metallicities. Also a low-Z star having the same ratio 
of surface velocity to critical velocity, w/wcrit (where Wcrit is the velocity for which the 
centrifugal force balances the gravitational force) as a solar-Z star has a higher surface 
rotation velocity due to its smaller radius (one quarter of Zq radius for a very low-Z 20 
Mq star). In the models presented below, the initial ratio w/wcrit is the same or slightly 
higher than for solar Z (see IHirschi 20071 for more details). This corresponds to initial 
surface velocities in the range of 600 — 800 kms~^. These fast initial rotation velocities 
are supported by chemical evolution models of Chiappini et al. (2006b) discussed in the 
next section. The mass loss prescriptions used in the Geneva stellar evolution code are 
described in detail in Meynet fc Maeder (2005). In particular, the mass loss rates depend 
on metallicity as M ^ [Z / ZqY'-^ , where Z is the mass fraction of heavy elements at the 
surface of the star. 

How do rotation induced processes vary with metallicity? The surface layers of massive 
stars usually accelerate due to internal transport of angular momentum from the core to 
the envelope. Since at low Z , stellar winds are weak, this angular momentum dredged 
up by meridional circulation remains in the star, and the star more easily reaches critical 
rotation. At the critical limit, matter can easily be launched into a keplerian disk which 
probably dissipates under the action of the strong radiation pressure of the star. 

The efficiency of meridional circulation (dominating the transport of angular momen- 
tum) decreases towards lower Z because the Gratton-Opik term of the vertical velocity 
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of the outer cell is proportional to 1/ p. On the other hand, shear mixing (dominating the 
mixing of chemical elements) is more efficient at low Z. Indeed, the star is more compact 
and therefore the gradients of angular velocity are larger and the mixing timescale (pro- 
portional to the square of the radius) is shorter. This leads to stronger internal mixing 
of chemical elements at low Z ( [Meynet &: Maeder 2002[ ). 

The history of convective zones (in particular the convective zones associated with 
shell H burning and core He burning) is strongly affected by rotation induced mixing (see 
IHirschi 2007p . The most important rotation induced mixing takes place while helium is 
burning inside a convective core. Primary carbon and oxygen are mixed outside of the 
convective core into the H-burning shell. Once the enrichment is strong enough, the H- 
burning shell is boosted (the CNO cycle depends strongly on the carbon and oxygen 
mixing at such low initial metallicities). The shell then becomes convective and leads 
to an important primary nitrogen production. In response to the shell boost, the core 
expands and the convective core mass decreases. At the end of He burning, the CO 
core is less massive than in the non-rotating model. Additional convective and rotational 
mixing brings the primary CNO to the surface of the star. This has consequences for the 
stellar yields. The yield of ^^O, being closely correlated with the mass of the CO core, is 
therefore reduced due to the strong mixing. At the same time the carbon yield is slightly 
increased. The relatively "low" oxygen yields and "high" carbon yields are produced 
over a large mass range at Z ~ 10~^ (jHirschi 2007^ . This is one possible explanation for 
the possible high [C/0] ratio observed in the most metal-poor halo stars (see Fig. 14 in 
[Spite et al. 2005| and IFabbian et al."2008|) and in DLAs (jPettini et al. 2008p . 

Models of metal- free stars including the effect of rotation (See contribution by Ekstrom 
and ' Ekstrom et al. 2008[) show that stars may lose up to 10 % of their initial mass due 
to the star rotating at its critical limit (also called break-up limit). The mass loss due to 
the star reaching the critical limit is non-negligible but not important enough to change 
drastically the fate of the metal-free stars. The situation is very different at very low 
but non-zero metaUicity ( [Meynet et al. 2006| iHirschi 2007) . The total mass of an 85 Mq 
model at Z = 10~* is shown in Fig. [T] (left) by the top solid line. This model, like 
metal-free models, loses around 5% of its initial mass when its surface reaches break- 
up velocities in the second part of the MS. At the end of core H burning, the core 
contracts and the envelope expands, thus decreasing the surface velocity and its ratio to 
the critical velocity. The mass loss rate becomes very low again until the star crosses the 
HR diagram and reaches the RSG stage. In the cooler part of the H-R diagram, the mass 
loss becomes very important. This is due to the dredge-up by the convective envelope of 
CNO elements to the surface increasing its overall metaUicity. The total metaUicity, Z, is 
used in this model (including CNO elements) for the metaUicity dependence of the mass 
loss. Therefore depending on how much CNO is brought up to the surface, the mass loss 
becomes very large again. The CNO brought to the surface comes from primary C and 
O produced in the He-burning region and from primary N produced in the H-burning 
one. 

The fate of rotating stars at very low Z is therefore probably the following: 

• M < 40 Mq-. Mass loss is insignificant and matter is only ejected into the ISM 
during the SN explosion (see contributions by Nomoto and Tominaga in this volume), 
which could be very energetic if fast rotation is still present in the core at the core 
collapse. 

• 40 Mq < Af < 60 Mq-. Mass loss (at critical rotation and in the RSG stage) removes 
10-20% of the initial mass of the star. The star probably dies as a black hole without a 
SN explosion and therefore the feedback into the ISM is only due to stellar winds, which 
are slow. 
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Figure 1. Left: Structure evolution diagram for a 85 Mq model with Vini = 800 kms ^ at 
Z = 10'^ Coloured areas correspond to convective zones along the Lagrangian mass coordinate 
as a function of the time left until the core collapse. The burning stage abbreviations are given 
below the time axis. The top solid line shows the total mass of the star. A strong mass loss 
during the RSG stage removes a large fraction of total mass of the star. Right: Composition in 
[X/Fe] of the stellar wind for the Z = 10"® models (solid lines). For HE1327-2326 {red stars), 
the best fit for the CNO elements is obtained by diluting the composition of the wind of the 40 
Mq model by a factor 600 (sec Hirschi 2007 for more details). 



• A/ > 60 Mq : A strong mass loss removes a significant amount of mass and the stars 
enter the WR phase. These stars therefore die as type Ib/c SNe and possibly as GRBs. 



3. Nucleosynthesis and galactic chemical evolution 

Rotation induced mixing leads to the production of primary nitrogen, ^^C and ^^Ne. 
In this section, we compare the chemical composition of our models with carbon-rich 
EMP stars and include our stellar yields in a galactic chemical evolution (GCE) model 
and compare the GCE model with observations of EMP stars. We also study the weak s 
process production that can be expected with the primary ^^Ne obtained in our models. 

3.1. The most metal-poor star known to date, HE1327-2326 
Significant mass loss in very low- 2" massive stars offers an interesting explanation for the 
strong enrichment in CNO elements of the most metal-poor stars observed in the halo 
of the galaxy (see Mevnet et al. 2006 [Hirschi 2007^ . The most metal-poor star known 



to date, HE1327-2326 (|Frebel et al. 2006P is characterised by very high N, C and O 
abundances, high Na, Mg and Al abundances, a weak s-process enrichment and depleted 
lithium. The star is not evolved so it has not had time to bring self-produced CNO 
elements to its surface and is most likely a subgiant. By using one or a few SNe and 



using a very large mass cut, Limongi et al. (2003) and Iwamoto et al. (2005) are able to 



reproduce the abundance of most elements. However they are not able to reproduce the 
nitrogen surface abundance of HE1327-2326 without rotational mixing. The abundance 
pattern observed at the surface of that star presents many similarities with the abundance 
pattern obtained in the winds of very metal-poor fast rotating massive star models. 
HE1327-2326 may therefore have formed from gas, which was mainly enriched by stellar 
winds of rotating very low metallicity stars. In this scenario, a first generation of stars 
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Figure 2. Chemical evolution model predictions of the N/0 and C/0 evolution, in the galactic 
halo, for different stellar evolution inputs. The solid curves show the predictions of a model 
without fast rotators at low metallicities. The dashed lines show the effect of including a popu- 
lation of fast rotators at low metallicities (for the data see |Chiappini et al. 2006al and references 
therein). 

(PopIII) pollutes the interstellar medium to very low metallicities ([Fe/H]^-6). Then a 
PopII.5 star (jHirschi 2005P like the 40 Mq model calculated here pollutes (mainly through 
its wind) the interstellar medium out of which HE1327-2326 forms. This would mean that 
HE1327-2326 is a third generation star. The CNO abundances are well reproduced, in 
particular that of nitrogen, which according to Frebel et al. (2006) is 0.9 dex higher in 
[X/Fe] than oxygen. This is shown in Fig.\l\{right) where the abundances of HE1327-2326 
are represented by the red stars and the best fit is obtained by diluting the composition 
of the wind of the 40 Mq model by a factor 600. When the SN contribution is added, 
the [X/Fe] ratio is usually lower for nitrogen than for oxygen. It is interesting to note 
that the very high CNO yields of the 40 Mq stars brings the total metallicity Z above 



the limit for low mass star formation obtained in Bromm & Loeb (2003) 



3.2. Primary nitrogen and ^^C 

The high N/0 plateau values observed at the surface of very metal poor halo stars require 
very efficient sources of primary nitrogen. Rotating massive stars can inject in a short 
timescale large amount of primary N. They are therefore very good candidates to explain 
the N/0 plateau observed at very low metallicity. According to the heuristic model of 
Chiappini et al. (2005) a primary nitrogen production of about 0.15 Mq per star is 
necessary. Using stellar yields calculations taking into account the effects of rotation 
at Z=10~® in a chemical evolution model for the galactic halo with infall and outflow, 
both high N/0 and C/0 ratios are obtained in the very metal-poor metallicity range in 
agreement with observations (see details in Chiappini et al. 2006a). This model is shown 
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in Fig. [2] (dashed magenta curve). In the same figm'e, a model computed without fast 
rotators (solid black curve) is also shown. Fast rotation enhances the nitrogen production 
by ^3 orders of magnitude. These results also offer a natural explanation for the large 
scatter observed in the N/0 abundance ratio of normal metal-poor halo stars: given the 
strong dependency of the nitrogen yields on the rotational velocity of the star, we expect 
a scatter in the N/0 ratio which could be the consequence of the distribution of the 
stellar rotational velocities as a function of metallicity. 

As explained above, the strong production of primary nitrogen is linked to a very active 
H-burning shell and therefore a smaller helium core. As a consequence, less carbon is 
turned into oxygen, producing high C/0 ratios (see Fig. 2). Although the abundance data 
for C/0 is still very uncertain, a C/0 upturn at low metallicities is strongly suggested 
by observations fsee lFabbian et al. 2008^ . Note that this upturn is now also observed in 
very metal poor DLA systems (jPettini et al. 2008)) . 

In addition, stellar models of fast rotators have a great impact on the evolution of the 
12(-i^13q ratio at very low metallicities (Chiappini et al. 2008). In this case, we predict 
that, if fast rotating massive stars were common phenomena in the early Universe, the 
primordial interstellar medium of galaxies with a star formation history similar to the 
one inferred for our galactic halo should have ^^C/^^C ratios between 30-300. Without 
fast rotators, the predicted ^^C/^^C ratios would be ~ 4500 at [Fe/H] = —3.5, increasing 
to - 31000 at around [Fe/H] = -5.0 (see Fig.2 in Chiappini et al. 2 008). Current data 
on EMP giant normal stars in the galactic halo ( [Spite et al. 2006] ) agree better with 
chemical evolution models including fast rotators. The expected difference in the ^^C/^^C 
ratios, after accounting for the effects of the first dredge-up, between our predictions 
with/without fast rotators is of the order of a factor of 2-3. However, larger differences 
(a factor of ^ 60 — 90) are expected for giants at [Fe/H]= —5 or turnoff stars already at 
[Fe/H]= —3.5. To test our predictions, challenging measurements of the ^^C/^^C in more 
extremely metal-poor giants and turnoff stars are required. 




Figure 3. Left: Structure evolution diagram (see description of Fig. [T| for a rotating 20 Mq 
model at Z = 10~® during H- and He-burning phases. Right: Chemical composition at the end 
of core He burning. Just above the core, one sees that the maximum abundance of ^'^Ne is 
around 1% in mass fraction and at the end of core He burning, around 0.5% is burnt in the core, 
providing plenty of neutrons for s process. 
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3.3. Primary ^^Ne and s process at low Z 

Models at Z = 10^^ show a production of primary ^^Ne during He burning. We also 
started calculating models at different Z to determine over which Z range the primary 
production of ^^Ne and also ^'^N is important. In Fig. [3l we show the properties of a 20 
Mq model aX Z — 10^^ up to the end of He burning. Around 0.5% (in mass fraction) 
of ^^Ne is burnt during core He burning and therefore leads to a significant neutron 
release. We calculated the s process by using the primary nitrogen of the rotating 20 Mq 
model aX Z = 10~^ inside a one-zone model (based on the Basel network and using an 
updated version of the reaclib reaction rate library) for s-process during He-burning. The 
first results are shown in Fig. 01 Large overproduction factors (100-1000) are obtained, 
however the process is not primary (see Pignatari et al. 2008 for more details). 
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Figure 4. Left: Ratio of the abundance at the end of He-burning to the initial abundance for 
a model including primary ^'^Ne. We obtain large overproduction factors between 100 and 1000 
in the mass range (A) 60 to 90. 



4. Conclusions and outlook 

The inclusion of the effects of rotation changes significantly the simple picture in 
which stellar evolution at low Z is just stellar evolution without mass loss. A strong 
mixing is induced between the helium and hydrogen burning layers leading to a significant 
production of primary ^*N, ^^C and ^^Ne. Rotating stellar models also predict a strong 
mass loss during the RSG stage for stars more massive than 60 Mq. These models 
predict the formation of WR and type Ib/c SNe down to almost Z=0. The chemical 
composition of the stellar winds is compatible with the CNO abundance observed in the 
most metal-poor star known to date, HE1327-2326. GCE models including the stellar 
yields of these rotating star models are able to better reproduce the early evolution of 
N/0, C/0 and ^^C/^^C in our galaxy. These new stellar evolution models predict a large 
neutron release during core He burning and we present here the corresponding s-process 
production during He burning. 

Large surveys of EMP stars (SEGUE, OZ surveys), of GRBs and SNe (Swift and 
GLAST satellites) and of massive stars (e. g. VLT FLAMES survey) are underway and 
will bring more information and constraints on the evolution of massive stars at low Z. 
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On the theoretical side, more models are necessary to fully understand and study the 
complex interplay between rotation, magnetic fields, mass loss and binary interactions 
at different metallicities. Large grids of models at low Z will have many applications, 
for example to study the evolution of massive stars and their feedback in high redshift 
objects like Lyman-break galaxies and damped Ly-alpha systems. 
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